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Abstract 

As well known, cumulant expansion is an alternative way to moment expansion to fully charac- 
terize probability distributions provided all the moments exist. If this is not the case, the so called 
escort mean values (or g— moments) have been proposed to characterize probability densities with 
divergent moments [C. Tsallis et al, J. Math. Phys 50, 043303 (2009)]. We introduce here a new 
mathematical object, namely the q—cumulants, which, in analogy to the cumulants, provide an 
alternative characterization to that of the (7— moments for the probability densities. To illustrate 
the technical details of the procedure, we apply this new scheme to further study a recently pro- 
posed family of scale-invariant discrete probabilistic models [A. Rodriguez et al, J. Stat. Mech. 
(2008) P09006; R. Hanel et al, Eur. Phys. J. B 72, 263 (2009)] having ^-Gaussians as hmiting 
probability distributions. 



PACS numbers: 05.20.-y,02.50.Cw,05.70.-a 



I. INTRODUCTION 



In classical thermodynamics for short-range-interacting systems we have two types of 
thermodynamical quantities, namely the extensive (e.g., total energy, total entropy, volume, 
etc) and the intensive (e.g., temperature, pressure, chemical potential, etc) ones. Within 
the formalism of Boltzmann-Gibbs statistical mechanics, these quantities typically emerge 
as successive moments of the A^-particle Hamiltonian. If, for whatever reason, exact calcu- 
lations are not tractable, a variety of procedures exist which involve truncations at some 
order. It is desirable that such truncations maintain the extensive or intensive nature of the 
quantities. Excepting the first moments in terms of the total Hamiltonian, which always 
are extensive, all the higher-order moments violate extensivity. There are however spe- 
cific combinations of these moments which preserve the extensivity. These are the so-called 
cumulants. For example, the total specific heat, which has to be extensive, appears as a 
second-order cumulant, i.e., a convenient combination of first- and second-order moments. 
Naturally, cumulant expansions are mathematically legitimate only when all moments are 
finite. A variety of physical systems exist for which this property is not verified. The pur- 
pose of the present paper is to develop a generalized form of cumulant expansion which 
overcomes this restriction. These generalized cumulants are in turn based on consistently 
generalized moments {escort moments or Q-moments), which we shall introduce later on. 
In Section II we introduce the g-cumulants; in Section III we review a special family of 
scale-invariant probability models; in Section IV we study the corresponding Q-moments; in 
Section V we address the associated (unnormalized) g-cumulants. We summarize our results 
in Section VI. 

Let us briefly remind the standard cumulant expansion. For a probability density function 
f{x) the cumulants kj = g^^\0), j = 1,2,..., are defined through the derivatives of the 
cumulant-generating function 



is the moment-generating function and fXj = (x^) = M^^\0), j = 0, 1,2, ... , are the mo- 
ments. Taking logarithms and Taylor expanding in Eq. (j2]) yields 



t^ 

g{t) = In M{t)=tKi + -K2 + ^K3 + -- - , 



(1) 



where 




(2) 
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InM(t) = ^121 



^3 

+ ^(/^3 - 3/il/i2 + 2/i?) 

+ ^(/X4 - 4/Xi/i3 + 12/iJ/X2 - 3/^2 - (3) 

5! 

+ 7T(Ai6 - lO/iq - 15yU2Ai4 + 30/i2 " 6yUiyU5 + 120/ii/i2yU3 + 30/i^/i4 

d! 

- - 120/i?/X3 + 360/i2/i4 - 120/i6) 



Comparing Eqs. ([3]) and ([T]) one gets ki = /ii, ^2 = /i2 — yU?, = yUs — 3/ii/i2 + 2/if 
which follow the general relation between moments and cumulants given by 



where the sum runs over solutions of the equation X]i=i ~ i' with Ui, i = 1,2, j being 
nonnegative integers, n = Yli=i ^'^d we have made use of the multinomial coefficient 
" ) = — — -.. Relation (H]) is a convenient re-writting of the relations given in 
l|, l2|. An alternative connection between cumulants and moments can be seen in |3|. 

As well known, the set of moments fully characterize a probability density function pro- 
vided they are all finite, the set of cumulants (linear combinations of the moments) being 
an alternative and, for some purposes (see later on), more convenient description. Once 
the set of moments are known, the probability distribution may be obtained via Fourier 
antitransforming the characteristic function <^(t) = M{it) = J-'[f]{t) = (ixe**^/(x). 

Notice that although K2 = {{x — (x))^) and K3 = ((x — {x))^), in general the cumulants 
do not coincide with the centered moments since ^4 7^ {{^ — (x))^) and so on. 

Let us review now an important property. We assume x = ^1+ ^2 + ■ ■ ■ + ^n, where {^1} 



are any N equal and independent random variables. We straightforwardly verify that 



{x) = iV(^i) 

{x') - {xf = - (6)'] (5) 

{x') - 3{x'){x) + 2{xf = - 3(e?)(6) + 2(^1)=^] 



In general 



n,iN) = NK,{l) (j = 1,2, 3,...; AT = 1,2, 3,...), (6) 



where the notation K,j{N) is self-explanatory. In other words, all cumulants are extensive 
in the thermodynamical sense. The main purpose of the present paper is to discuss what 
happens with this property in the presence of strong correlations such as those that are 
typical within g— statistics [3], a current generalization of Boltzmann-Gibbs statistical me- 
chanics (recovered as the q = 1 particular instance). This generalization, sometimes referred 
;o as nonextensive statistical mechanics, has received a wide variety of physical applications 



19|. 



II. g-CUMULANTS 



Within the frame of g— statistics, a generalization of the concept of cumulant is in order. 
A first natural attempt for such generalization consists in just replacing the logarithm and 
exponential functions in ([3]) by the g— logarithm and g— exponential functions (see Appendix 
for definitions and corresponding expansions). Their respective Taylor expansions yields 
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, / txv . 

t 

= 

+ ^^((2^ - l)/^3 - 3g/ii/i2 + (g + 1)^?) 

+ ^g((2g - l)(3g - 2)fi, - 4g(2g - l)/ii/i3 + 6g(g + l)/i?/^2 - SgV' + 2)/^f) 

^5 

+ -g((2g - l)(3g - 2)(4g - 3)^5 - 5g(2g - l)(3g - 2)^i/i4 - 10g2(2g - l)/i2/i3 
5! 

+ 10(g + l)g(2g - 1)^2^3 - Wq{q + l)(g + 2)/i?/i2 + 15g2(g + 
+ (g + l)(g + 2)(g + 3)^?) 

: (7) 



where we have defined the q—cumulant of order j, Kj^q, which is related with the moments 
through 



where e{q,n) and l{q,n) are the coefficients of the Taylor expansion of the g— exponential 
and g— logarithm functions respectively (see Appendix). 

Notice that the moments {fij} in ([7]) are exactly the same introduced in Eqs. ([2]) and Q. 
In other words, {fij} do not depend on q. Notice also that Ki g = /ii, Vg (the g— cumulant 
of first order coincide with the first moment for any value of q), and Kj^i = Kj Vj (the 
g— cumulants reduce to the standard cumulants for g = 1). 

Although interesting in principle, we will show that the g— generalization which is useful 
is not the above one, but the one which follows now. We will define the unnormalized 
g— cumulants, k'^ ^, through 
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^3 ^4 ^5 

_t , ' 

+ ^^((2? - l)/^3,3g-2 - 3g/i2,2g-l/^'l,g + + 

+ ^^((29 - l)(3g - 2)/i'4^4^_3 - 4g(2g - l)/i3,3q-2/^'l,q + 6g(g + l)/^2,2q-l/^'L 

-3gV2,2,-i-(? + l)(? + 2)/it,) (9) 

+ ^-q{{2q - l)(3g - 2)(4g - 3)/i'5,5,_4 - 5g(2g - l)(3g - 2)^^',^,^_,^^',^^ 

- 10g2(2g - l)M3,3g_2/^2,2g-i + 10(g + l)g(2g - l)/i3,3,_2/^1,q 

- 10g(g + l)(g + 2)^'2,2g-i/^t, + 15g2(g + 
+ (g+l)(g + 2)(g + 3)/xt,) 

where we shall call Mq{t) = {e'q^'^^^'' ) the Q-moments generating function. Expression ([9]) 
makes use of the so called unnormalized escort moments [20], ^'jq = J dx f (x)'^ , related 
in turn to the standard escort moments (or Q— moments), fij^g = l^'j^q/ J dxf{x)^, i. e., 
x^ averaged over the escort distribution /q(x) = /(x)*5/7q, with 'Jq = J dx f{x)^. In the 
limit Q — 1, the escort distribution as well as the Q— moments, and the unnormalized 
Q— moments, tends to the standard ones. 

Notice from Eq. ([9]) that, for the first unnormalized g— cumulant we must use the 
Q— moment with Q = g (in fact the unnormalized g— cumulant coincides with the un- 
normalized escort moment for j = 1: k,[ ^ = fi[ g), for the second unnormalized g— cumulant, 
Q— moments with two different values, namely Q = q and Q = 2g — 1, must be used, and sim- 
ilarly for higher-order unnormalized g-cumulants. The relation between the unnormalized 
g— cumulants and the unnormalized Q— moments is thus given by 

"- = ^' ^ ^^^(^'"\n,n2 ^ I (10) 

ELi*"«=i J' i=i \ J 

For g = 1, the unnormalized escort moments appearing in Eqs. ([9]) and ( ITOj) : /i'jg, with 

Q = ^(9 ~ 1) + 1; ^ = 1) 2, 3, ... , reduce to the normalized escort moments which in turn 
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coincide with the ordinary moments: /i'j ^ = /ij^i = fii The same happens with the 
unnormahzed g— cumulants: i^'n = = Wi, so Eqs. ([8]) and ffTOl) coincide. 

For the general q ^ 1 shown in Ref. j2o|, the family of unnormalized es- 



cort moments /i! 



= 1,2,3,..., provided they are all finite, may serve, together 
with their acompanying denominators, 7j(q-i)+i = / dx f{xy^'^~^^~^^ , in the expresion of the 
corresponding escort moments, as an alternative description for the quite frequent cases 
of probability density functions whose moments are not defined (e.g. moments of order 
j ^ 2 for q—Gaussians with q ^ 'b/?) 2l|]). In analogy with the Q = 1 case, once the 
Q— moments generating function is known, the probability distribution may be obtained 
via g— Fourier antitransforming the q- characteristic function ^q{t) = Mq{it) = J^q[f]{t) = 
dx eq^^^^^''__f (x) , (/(x) ^ 0, g ^ 1), where we have made use of the so called q-Fourier 
transform 22|,|23|, which is a (nonlinear) g— generalization of the standard Fourier transform. 



III. SCALE INVARIANT TRIANGLES 

We will apply the above formalism to further study a family of scale- invariant probabilistic 
models characterized by a real number u > 0, first introduced in [2J] and later generalized in 



251]. The model consists in a set of equal, long-range- correlated binary random variables 
{xi}, which take values and 1. 

For the 2^ elementary events of the sample space of the N variables system there is a 



set of only (A^ + 1) different probabilities given by 25| 



H _ B{n + u,N -n + u) 



for n = 0,1, . . . , N, where B(x, y) stands for the Beta function. The {?")v „} satisfy the so 
called Leibnitz triangle rule 24] 

I „M _ „H /-ION 

' N,n "T ' N,n+l ~ ' N-l,n V^-^J 

which is a fingerprint of the scale invariant character of the system, which states that the 
marginal probability distribution of the N system coincides with the joint probability dis- 
tribution of the corresponding (A^ — 1)— subsystem. In other words, in its continuous ver- 
sion, we have J Pn{xi,X2, ■ ■ ■ ,XN)dxN = Pn-i{xi,X2, ■ ■ ■ ,xn-i) (although similar in form, 



this conditions are not the Kolmogorov consistency conditions of a stochastic process 



see, for instance 26| — . Indeed, the Kolmogorov conditions refer to a system of ele- 
ments whereas here we impose nontrivial conditions connecting the joint probabilities of a 
(A^ — 1)— subsystem with the marginal ones of a A^— system). 

For a fixed value of z/, probabilities f llip may be displayed as a function of in a 
symmetric (i.e., r^|^ = T^^jy-n) triangle (thus the name scale invariant triangles). As an 
example, the = | triangle reads 



{N = 0) 
iN = l) 
{N = 2) 
(AT = 3) 
(AT = 4) 
(AT = 5) 



1 1 

3 ^ 1 ^ 3 

1 ' 1 ' 1 ' 1 

35 1« 5 1« 3 1« 5 1« 35 



63 7 3 7 1 ^2*^ 63 



256 256 256 256 256 256 



Due to the properties of the Beta function the coefficients of the triangle r^p^ are rational 
numbers for any rational value of u. 



Special values of u are z/ = 1, for which the original Leibnitz triangle 27| is obtained 
[l/r^Nn = {N + 1)(^)], and the limiting case r^^^ = lim;^_>oo = 1/2^, for which the 
system becomes uncorrected 



2J]. In the following, we shall refer to this limiting 
the Leibnitz and the de Moivre-Laplace (or Boltzmann-Gibbs) limits respectively. 

The probability distribution for the discrete variable z = xi + X2 + ■ — h xn, which takes 
values 0, 1, . . . , A^, is given by 

pZ^Piz = N-n)=(^^^y;;l (13) 

with X^^o^'iVn ~ to degeneracy imposed by the identical character of de A^ 

binary variables. By properly centering and scaling with a map n — )■ y{n), the distribution 
f|T3|) is transformed to a new one, V^'^^y), which tends to a g— Gaussian 

P('^)(2/)^ (14) 



for — > oo, with Ag such J A^e'^^dy = 1, and a value of g = qum which depends on u and 
the chosen change of variable. The variable change, the intermediate distribution, and the 
value of giim are given in Table [11 The Leibnitz case [u = 1 hence q — t- — oo) corresponds to a 
uniform distribution while for the Boltzmann-Gibs case {u = oo hence q = 1) the standard 
Gaussian is obtained. 



map 



y = 2V^7^ ( ^ - i 



n 



n) 




4(n(iV 



iV2 



1 

V + - 

2 



< 1; {j^> 1) 



2u + 3 
2v + l 



> 1; {v> 0) 



TABLE L The two different changes of variables 



giim < 1 
and giim 



24 



corresponding respectively to all real values 



25] and all real values 1 < giim < 3 25], and the associated the corresponding V^'^\x) 



Thus, depending on the chosen discretization (namely the Table I top or bottom dis- 
cretization), for each value of z/ ^ 1 there exist two conjugated g— Gaussians, one of them 
with 1 > q = qu E (— oo, 1) (with compact support \x\ < ^^j-^ ) and the other with 
1 < g = G (1,5/3] (the support being the whole real axis). In the case of u being a 
positive integer, the different values of g are shown in Table HTl For u G (0, 1), however, only 
one g— Gaussian with 1 < q = q^ E (5/3, 3) is obtained 25|. Figured] plots the two different 
giim values as a function of z/. 

Relations q{v) given in the third column of Table [I] can be inverted so as to express the 



u 


0+ 1/2 


1 2 3 4 •• 


• oo 


qu 




-oo 1/2 2/3 •• 


• 1- 




3- 2 


5/3 7/5 9/7 11/9 •• 


• 1+ 



TABLE II. Values of g as a function of for the two families of conjugated (7— Gaussians. 
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family of triangles function of q |25| : 



) 



g < 1 



q = l 



(15) 



B (n + 



/V - n + 



2{g-l) 



) 



1 < g < 3 



B 



( 3-g 3-1? 



) 



V 



\2{q-l)^ 2{q-l) 



which allows as to fix a priori the value of q and then get the triangle generating the 
corresponding g— Gaussian. 

Figure E] shows the probability distribution V^'^\x) ioi a N = 1000 system and two dif- 
ferent values of the parameter u, namely z/ = 2, for which gj, = < 1 and qu = 7 /5 > 1, and 
z/ = 1/2, for which there exist one only value q^, = 2. The corresponding g— Gaussians (not 
shown) graphically overlap with the shown distributions. To further check Eq. f|T^ . Fig. [3] 
shows logg^(p('')(x)/'p('')(0)) versus (xp('')(0))2 — which, in the — )■ oo limit, coincides 
with a straight line for a g— Gaussian — for u = 1/2, q^ = 7/5 (left) and u = 2, q,^ = 2 
(right) for = 200, 500 and 1000. As can be seen in the figure, points far from the origin 
(corresponding to small values n = 0, 1, . . . , and their symmetric n = A^ — 1,A^ — 2,...) devi- 
ate from the fitting to a straight line. Nevertheless this is a finite size effect since increasing 
A^ has the effect of bringing them back to the line. For values of u and discretizations corre- 
sponding to compact support g— Gaussians (not shown) this finite size effect is not present 
and the fitting is improved. 

IV. SCALING OF THE Q-MOMENTS 

Our aim in this section is to study the scaling with the system size of the Q— moments 
and g— cumulants of the family of triangles (ITT]) or, equivalently (fTSjl . 
For our model, the escort moments of interest read 




(16) 



n=0 



10 



where the numerator of the last fraction stands for the corresponding unnormahzed Q— moment, 
the denominator being the normahzation factor of the escort distribution. 

Straightforward calculations based on the symmetric character of triangles fill I) allows us 
to exactly calculate the first-order escort moments as 

M:i=4:i=f, vg, vz.. (17) 

The higher-order escort moments for the aforementioned Boltzmann-Gibbs and Leibnitz 
limiting cases are addressed next. 



A. Boltzmann-Gibbs limit 



To begin with, we shall study the standard {Q = 1) moments for the Boltzmann-Gibbs 
{u — 7- oo) limit. In this case the probability distribution f lT3|) reduces to the binomial 
distribution pj^j = 2^(^) (which tends to the Gaussian distribution in the thermody- 
namic limit — > oo). Taylor expanding its moment generating function M{t) = (e*^) = 
2^En=o = ((1 + e*)/2)^, the standard moments n^^^ = M^^\0) of the binomial 

distribution are easily obtained: 



(oo 



(oo 



(oo 

(oo 
f^4A 



(oo 
/^5,1 



N{N + 1) 
4 

N^{N + 3) 
8 

N{N^ + 6N^ + 3N -2) 
16 

N{N^ + lON^ + 15A^2 - ION) 
32 



(18) 



Though not shown here, it is worthy noticing the unexpected result that the scaling law 
( ITSl) holds for the escort moments of the binomial distribution for any value of Q, that is 
for the binomial distribution one gets 



(oo) 



1 

2^ 



(19) 



On the other hand, the standard {q = 1) cumulants k^'^'* of the binomial distribution 
can be obtained making use of Eq. (jl]) or, alternatively, as the derivatives of the cumulant 
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generating function g{t) = A^ln((l + e*)/2): 



K 



(~) 
1.1 



N 



K 



2,1 



N 
1 



— , K, 



(~) 
3,1 



(oo) 
M,l 



"5,1 



(20) 



where all the odd cumulants vanish from the third on. It has to be noticed that, for 
—7- oo, while the moments follow the rule /ij'^'' oc A^-', the cumulants remain extensive, 
that is, oc A^. In fact, we easily verify that the results in Eq. f l20|) satisfy Eq. ([6]). 



B. Leibnitz limit 



In the specially simple Leibnitz case (z/ = 1), the probability distribution ( 1131) reduces to 
the uniform one = which once introduced in Eq. (lT6l) yields /i^^q = X]n=o^"' 
for all Q, so the escort moments are independent of the value of Q, thus coinciding with 
the standard {Q = 1) moments, which may again be calculated as the derivatives of the 
corresponding moment generating function M{t) = (e*^') = jy+iXl^Lo^*" ~ ]v+T^~~rS~~"' 
M(0) = 1: 



Pl,Q — 


N 
~2 






u^') - 




N 






X 










ATS 

+ 3 — + — 
10 30 


N 
" 30 


f^5,Q — 


~~6~ 


iV3 

+ - + T2- 


' 12 



(21) 



In turn, the standard [q = 1) cumulants k^^] of the discrete uniform distribution can be 
obtained through Eq. (jlj) or from the cumulant generating function g{t) = \n ' ^ i-e^^+i)* 



, N+l l-e* 

12 



g{0) = 0: 



K 



(1) 
1,1 



(1) _N\N 
''^'i "12+6 



K 



K 



K 



(1) 

3,1 

(1) 
4,1 

(1) 

5,1 









120 30 



(22) 



iV2 

20 



N 
30 



where again the odd cumulants of order greater than one vanish. Scahng of cumulants ( 122|) 
may be expressed as nfl ~ Bj being the Bernoulh numbers {Bi = 1/2, B2 = 1/6, 

B3 = 0, 54 = -1/30, B5 = 0,Be = 1/42, . . . )• 

C. General case 

We shall now allow u to take on any positive value in (fT6l) . In principle, Q may take on any 
real value. However, as we will be interested in the escort moments involved in the expression 
of the unnormalized g— cumulants fllOp . we will first focus on the case Q = j'{q — 1) + 1 (the 
extra index j' has been introduced in order to allow for a more general description though 
the only escort moments appearing in Eq. ffTOl) are those with j' = j). Furthermore, for the 
escort moment of the u triangle we will take either q = ot q = as given in Table [H 

Figure m shows the scaling of the escort moments of orders j = 2, 3 and 4, for 

u = 2 and 3 and different values of j'. The following trend is observed 

~ Vz., Vj' (23) 

so the exponent of the scaling of the escort moments depends only on the order j of the 
moment not on j', neither on q^, which considerably simplifies the discussion. From now on 
we will restrict ourselves to the case j' = j (and correspondingly denote Aj^g^ = Ajj/=j^q^ in 
fl23|) ). Not shown results reveal that Eq. fl23l) is also valid for Q = j'{q — 1) + 1 with q = q^. 

The prefactor Aj g^ in the saling law (123|) has already been analitically obtained in the 
Leibnitz (g = —00) case as ^j,-oo = as well as in the in the opposite limit, the Boltzmann- 
Gibbs (g = 1) case as Aj^i = ^ (see Eqs. (l2Tll and (|T8l) respectively). A transition between 

13 



both behaviours is expected for the hmit g — )■ 1^. In order to numerically obtain Aj^q^, 
Fig. [5^ shows the ratio f^^j'j(^g for = 50 as a function of for different values of 

j. It is observed that Aj^g^^ ~ jj. Figure shows again Ajg^, now approximated as the 
ratio fJ'^^j(^g^_i^^i/N^ for N = 500 as a function of u. It is clearly observed the expected 
asymptotic limit 

lim Aj^g^ = y (24) 

A further generalization for the scaling of the Q— moments must be done. Not shown 
calculations reveal that relation (!23|) still holds when Q = j{q — 1) + 1 and q takes on any 
real value, not necessarily the values indicated in Table [H (which in turn implies that Q may 
take on any real value). Thus, it can be stated that 

which means that the whole family of triangles ffTTj) . or equivalently f|T5|) share the same 
scaling behaviour for their escort moments /Wj'q- 



V. SCALING OF THE UNNORMALIZED g-CUMULANTS 

We shall now turn to the study of the unnormalized g— cumulants k'-^ given in Eq. ffTOj) . 
Insomuch the unnormalized g— cumulant of order j, n'- g, is a linear combination of products 
of unnormalized Q— moments, which is homogeneous in the order j of the g— cumulant (see 
restriction X]i=i^'^i — J sum defining the unnormalized g— cumulant in Eq. ( fTOl) ). 

the scaling of n'jg will coincide with the scaling of the unnormalized Q— moment fi'jQ with 
Q=j{q-1) + 1. 

Fi gure|6]shows the asymptotic behaviour of the unnormalized escort moments f^' j j(q-i)+i 
for the z/ = 3 (gs = 1/2, gs = 9/7) triangle for g = 3/2 (left) 2 (center) and 5/2 (right). The 
following scaling relation is fulfilled 

Thus, the independence of the value of g characterizing the scaling exponent of the 
escort moments (see Eq. ( |25l) ) no longer holds for the unnormalized escort moments. Some 
comments must be made on the validity of scaling law (!26|) . It is valid for any positive g 
for ^ 1 (though it may fail for g close to and values of j > 1). In the case u < 1, 

14 



Eq. (1261) holds only for q ^ 1, while the exponent scaling depends on the value of z/ < 1 
in a complicated fashion for g > 1. It must be said, however, that relation f l26|) is only 
approximated when q approaches 2. For g = 2, in the transition from the increasing to the 
decreasing trend of f^'^^g with A^, a slight increase of the unnormalized Q— moments with 
is observed (see Fig. E]) instead of the expected constant value. We will discuss this point in 
detail later when studying the scaling of the unnormalized g— cumulants. Finally, for q = 1 
the unnormalized moments reduce to the standard ones and relations ( 125|) and fl26|) coincide. 

By comparing Eqs. (fT6|) . ( 125|) and ( |26|) . it is readily deduced the scaling law for the 
normalizing coefficients of the escort distributions 7q as 

7]ti)+i - f^AT^-^-^) (27) 

within the same range of parameters as stated above. For the q = 1 case, Eq. fl27j) is 
trivially fulfilled, representing the normalization of the family of probability distributions 

As already mentioned, an analogous relation to (l26l) holds for the unnormalized g— cumulants 



12 - C,,,,.N^^'-'^^ (28) 



with the same restrictions for the values of g and u. Special cases are g = 0, for which all the 
unnormalized 0— cumulants vanish except for the first one (see Eq. (jO])), which coincides 
with the corresponding unnormalized escort moment n'^^l = fi''^^ oc A^^, and g = 1, for 
which the standard cumulants are recovered as well as an analogous scaling law as that for 
the Leibnitz case f l22|) for any value of z/, that is, k^j^I oc A^-' with vanishing odd cumulants 
from the third on. 

Nevertheless, relations (l26l) to ( !28l) are no longer valid in the z/ — )■ oo limit, that is, for 
the binomial distribution. It may be shown that in this case the scaling for A^ — > oo is 

oc oc Ar^'(3-9)/2. (oo) ATj(l-q)/2 (2n\ 



which coincide with relations (1261) to ( 1281) for g = 1. The very slow and nonuniform character 
of the convergence in limi,_!.oo t^|^ = 1/2^ and the numerical difficulties in dealing with 
coefficients r^^„ for large u and A^ make it subtle to study in more detail the transition of 
the scaling exponents from the z/ < oo to the Boltzmann-Gibbs case. Figure [7] shows the 
scaling of the first unnormalized 3/2— cumulant for triangles with u = 10, 20 and 30, together 
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with the Boltzmann-Gibbs case. The crossover from the scahng exponent 1/2 given by (128|) 
for < oo and exponent 3/4 indicated in fl29|) for z/ = oo is clearly observed. Finally, Eq. f l29|) 
is valid for g ^ 1 for the unnormalized escort moments and their corresponding 

normalizing coefficients 7-?° "'in,,. However, for the unnormalized o— cumulants k'^-^^ its 
validity restricts to g > 1. The case q = 1 must be excluded since the 1— cumulants 
(normalized or not) coincide with the standard ones whose scaling k^j^^ oc has been 
already obtained in fl20|) . 

Let us connect now the values of u and q in Eq. (!28|) . i. e., let us study, for the z/ triangle 
f fTTj) . the g— cumulants with g being the corresponding limiting value of the g— Gaussian 
( IT^ . As stated in Table HI there is only one limiting value, gum = q^ for u E (0, 1), but two 
different options, namely giim = qu and giim = g,y, for z/ ^ 1. We shall denote with g(,^) any 
of them and calculate the unnormalized g— cumulant k'^-^} . 

Following the comments made for the validity of Eqs. fl26|) to (l28l) . we may substitute 
g by g(jy) in Eq. fl28|) if u ^ 1 since it is valid for any value of g. The u < 1 case must be 
studied separately since for this case q(^^) = q^ > 1 and, as already mentioned, Eq. (!28|) no 
longer holds. Nevertheless, provided 1/2 < z/ < 1, Eq. ( l28i) still holds but now only for the 
specific value g = g(,^). On the contrary, for z/ < 1/2 (hence g;^ > 2), the scaling exponent 
reverts its sign. Sumarizing, the following scaling law is obtained 

-It) - (30) 

Some comments on the validity of relation ( l30l) must be made. As shown in Fig. |H1 for 

(-) {-) 

gi = 2, a logarithmic increase with A^ of kW = A^'i 2 is obtained instead of the expected 

constant value predicted by (1301) . which therefore is no longer valid in this case and has to be 

(-) 

replaced by a logarithmic law ~ A+Bln N whose coefficients are found to be A = 0.315 
and B = 0.105. This logarithmic correction influences numerically the scaling relation ( 130|) 
for g(,^) close to 2, as can be seen in Fig. [81 Finally, for the g < case, the scaling law still 
works (150]) for j = I. Figure [H] shows the actual value of the scaling exponent Q;i,g(^j = |2— g(,y) | 
of the flrst order unnormalized g— cumulant k,'['^1 = ul'f}, ~ Ci o, , A^^^'''^") . A small 
deviation from the scaling law f l30|) is observed for qi^y) > 2. 

Concerning the proportionality coefficient in (!30|) . it must be noticed that 
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so the dominant term in the scahng law (l30i) vanishes for j > 1. Subdominant terms are 
also negligible, but in the limit z/ — t- oo a dominant linear term arises so one recovers the 
extensive behaviour k^°^'^ ~ advanced in (I2U]) . 

Left panel of Fig. [TO] shows the order 1 cumulants in the Boltzmann-Gibbs limit 
1/ — )■ oo. The predicted trend lim,^_>.oo I'^i'^q^^^ = i^iT^ = y is clearly observed. Right panel 
of Fig. [TD] shows the limit flHTl) for j = 1 and 2. (Notice that Figs. [7] and [10] differ in the 
fact that for Fig. [7] we take the u oo limit for a g— cumulant with a fixed q = 3/2 value 
whereas in Fig. [T^we take the same z/ — )■ oo limit but for the corresponding g(i,)— cumulant, 
so g(^) 1). 



VI. CONCLUSION 

We have studied in detail the scaling with the system size of the Q— moments as well as 
the g— cumulants, — normalized or not — , of the family of scale invariant triangles introduced 
in Section lllli We summarize the scaling laws that we have established in Table IIIIl The 
prefactors have been indicated whenever known. For the range of validity of expressions in 
the last column of the Table, see Section [V] The scalings with N of the u < oo column 
reduce to the z/ — )■ oo column for g = 1 except for the unnormalized g— cumulants for j > 1. 
This unexpected feature (see Fig. 7), i.e., the fact that the g — )■ 1 and N ^ oo limits do not 
commute, constitutes an interesting result of the present paper. 

We believe that the present generalization of the cumulant expansion could be of interest 



for analytically discussing long-range- interacting systems (e.g., 28|, |29|]), for which moments 
above a given order typically diverge. 



Appendix 

For g G M, g 7^ 1, the g— exponential function and its inverse, the g— logarithm function 
are defined as 

e^^[l + (l-g)x]T^ (A.l) 

- 1 , , , 

\ngx = — (A.2) 
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TABLE III. Unnormalized Q— moments, normalization factors, ordinary Q^moments and unnor- 
malized (7— cumulants for de Boltzmann-Gibbs case (third column) and the family of triangles 
(fourth column) for j = 1, 2 and the general case. In all cases, the symbols ~ and oc are to be 
understood in the N ^ 00 limit. 



where the standard logarithm and exponential functions are recovered in the limit o — )• 1. 
Their respective Taylor expansions are given by (see, for instance, the last of Refs. |^) 



n=0 



e{q,n) 



1; 



n = 0,l 



g(2g - 1) ■ ■ ■ ((n - l)g + 2); n>l 



\n— 1 



n=l 



-l{q,n)x''; l{q,n) 



1; n = 1 

g(g + l)---(g + n-2); n> 1 



(A.3) 
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FIG. 1. Limiting values Qi, and q,^ as a function of u. Notice the u^q) is a single- valued function, 
whereas q{i') is not. 




-4-2 2 

X 



FIG. 2. Probability distribution V'^'^^x) for = 2, g,, = (solid line), = 2, q„ = 7/5 (dotted 
line), and v = 1/2, q^ = 2 (dashed line). = 1000 in all cases. 
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FIG. 3. logg^{r^''\x)/r^''\0)) versus {xr^''\0)f for z/ = 1/2, = 7/5 (left) and u = 2, q^, = 2 
(right) for N = 200 (triangles), 500 (crosses) and 1000 (circles). The dotted straight line is obtained 
by linear fitting after dropping the last four points. 
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FIG. 4. Scaling of the escort moments ()16p with Q = j' {qu — 1) + 1 for different values of j', for 
V = 2 (up), u = 2i (down) and j = 2 (left), j = 3 (center) and j = 4: (right). The behavior for q = Qu 
is totally analogous. For the not shown j = I case all the curves collapse in one: _i)+i = T 

since relation (fT7|) is satisfied. 
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FIG. 5. a) Ratio fJ-jjf^g with i/ = 50 as a function of N and j = 1,2,3,4 and 5 from top 

to bottom, b) Aj^q^, calculated as the ratio /^^'^^•(g with = 500, as a function of u and 

the same values of j as in a) . 
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FIG. 6. Scahng of the unnormahzed escort moments f^'^j^j(q~i)+i triangle = 3 (q^ = 1/2, 

qs = 9/7) for q = 3/2 (left) 2 (center) and 5/2 (right) for j = 1 to 5. The scaling law ([26]) is 
obtained. 
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FIG. 7. Scaling of 

'^'13/2 — '^'i 3/2 ^ — ^^'^ '^'^^ scaling exponent is 1/2 for 

< 00 and 3/4 for 1/ = 00. 
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FIG. 8. Log-log plots of k'i ^ = fJ- iq^ ^ versus AT for = Qv = 1-6 to = 1.9 (left) and Qi, = 2.1 
to Qi, = 2.4 (right). The center panel shows the qi = 2 case and its fitting to a A + 5 In A?^ curve 
with A = 0.315 and B = 0.105. Both curves virtually overlap. 
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FIG. 9. Scaling exponent of the first order unnormalized g— cumulant ^'i*^,]^ ^ = [^'\\^ j ~ 

^1.9(1.)-^"^ • ^ small deviation from the scaling law ai,q(^,) = |2 — given in Eq. (pO]) is 
observed for g^^) > 2. Exponents have been obtained via linear regresion of the ln(K'^'^g^ ^) versus 
In curves for N from 500 to 1000. The \q — 2| curve is also plotted as a guide to the eye. 
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FIG. 10. (Left) Log-log plot of n i ^ = A* i j versus for conjugated pairs of triangles = 2, 3, 
10 and 20. The limiting curve = N/2 is plotted as a guide to the eye. (Right) Proportionality 
coefficient Cj^q^^^ in ([30|) for j = I and 2 as a function of q(^,yy The expected limit (f3T]) is obtained. 



27 



